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[\\  This  paper  describes  a  method  for  estimating  spatially  variable  bottom  roughness 
lengths  (A>)  in  friction-dominated  coastal  regions  where  dense  measurements  of  the 
significant  wave  height  are  available.  The  method  utilizes  a  numerical  wave  model  to 
calculate  wavefields.  The  model -predicted  significant  wave  height  is  compared  to  a 
control  simulation  with  a  known  kb  field,  which  is  a  proxy  for  measured  wave  heights.  The 
error  is  used  in  combination  with  an  influence  matrix  to  successively  correct  the 
bottom  roughness  field.  This  predictor-corrector  calculation  is  completed  in  a  series  of 
analysis  cycles.  The  method  is  demonstrated  in  an  idealized  basin  with  different  kh 
distributions.  The  test  cases  simulate  swell  propagating  over  a  sloping  beach.  The  original 
kb  fields  are  recovered  in  a  reasonable  number  of  analysis  cycles  but  the  method  is 
limited  by  the  influence  of  bottom  friction  on  the  wave  height.  The  inversion  is  shown  to 
be  robust  in  the  presence  of  errors  in  the  measured  wavcfteld  as  well  as  random 
bathymetry  errors.  However,  the  inversion  fails  if  bathymetry  errors  are  large  and/or 
systematic  because  the  friction  error  is  not  substantially  greater  than  the  error  from 
bathymetry,  which  is  also  a  key  parameter  for  calculating  the  wave  height.  Thus  it  is 
important  to  select  parameters  and  variables  that  have  well-defined  dependencies  in  the 
numerical  wave  model  for  this  procedure  to  be  effective. 
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distributions  u>ing  a  numerical  wave  model,./.  Geophys.  Res. ,  112.  C08008.  doi :  10/1029 '2005 JC0G3 309. 


L  Introduction 

[:]  Jt  is  not  possible  io  completely  observe  the  global 
ocean  but  remote  sensing  has  made  it  easier  to  measure  sea 
surface  properties  over  large  regions.  The  observation 
problem  is  even  greater  for  ihe  seafloor.  Although  bottom 
friction  is  not  critical  to  understanding  the  dynamics  of 
ocean  currents  in  deep  water,  it  can  be  important  for  both 
currents  and  waves  in  ihe  shallow  water  over  the  continental 
shelf.  When  the  water  depth  is  less  than  the  deep-water 
wavelength  of  wind-generated  surface  waves,  bottom  friction 
can  produce  observable  changes  in  the  surface  wave  proper- 
lies  (most  obviously,  dissipation  of  wave  energy  as  waves 
propagate)  as  demonstrated  in  recent  work  [Sheremet  and 
Same,  2003;  Ardhuin  et  al 2003:  Kaihatu  and  Sheremet, 
2004].  These  changes  can  be  observed  remotely,  as  can 
other  ocean  surface  properties  like  wave  energy  spectra  and 
temperature  However,  when  only  surface  properties  are 
available,  it  is  problematic  to  determine  the  causes  of  the 
observed  changes.  Because  of  these  difficulties,  ocean 
scientists  have  begun  to  utilize  inverse  techniques,  which 
can  improve  our  knowledge  of  physical  processes  from 
observations. 
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[3]  One  kind  of  inverse  method  is  data  assimilation, 
which  combines  model  physics  with  observations  to  pro¬ 
vide  a  better  picture  of  the  ocean  than  can  be  deduced  from 
either  alone  [Anderson  el  al,  1996],  Data  assimilation 
techniques  range  from  nudging  numerical  models  with 
observations  to  direct  assimilation  of  observations  using 
variational  approaches  |Lc  Dime!  and  Tala  grand,  1986; 
Bennett,  1992;  Bertino  et  al,  2003],  Data  assimilation  has 
been  used  lo  improve  numerical  wave  forecasts  by  nudging 
wave  models  w  ith  observations  from  wave  buoys  and 
remote  sensing  [O  Rciflv  and  Guza,  1998;  Hnhhuijsen 
et  ai%  1997;  Bidlot  and  Holt,  1999;  Greenslade,  2001],  In 
addition  to  its  use  in  improving  wave  forecasting,  assimi¬ 
lation  has  also  proven  useful  in  estimating  water  depth  using 
inverse  techniques  [Dahymple  et  al ,  1998;  Grills,  1998; 
Wackerman  et  at. ,  1998],  and  iterative  approaches  are 
necessary  for  depth  inversions  based  on  numerical  wave 
modeling  [Kennedy  et  al,  2000]. 

[4]  Data  assimilation  requires  that  some  assumptions  be 
made  wiih  respect  lo  the  relationships  between  the  data  and 
model  parameters  being  investigated.  Numerical  wave 
models  lira!  include  physical  phenomena  such  as  shoaling, 
dissipation,  bottom  friction  and  refraction  can  permit  a 
comprehensive  examination  of  complex  processes.  However, 
in  order  to  use  a  numerical  wave  model  in  an  inverse 
solution,  it  is  necessary  to  first  identify  the  key  model 
parameters.  Key  model  parameters  are  those  upon  which 
all  other  parameters  have  a  high  dependency.  They  can  be 
identified  by  sensitivity  analyses  using  a  numerical  model 
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[Sun  et  at.,  20(1 1 ;  WWxse  ami  Fesen  2003)  but  it  is  important 
to  test  for  consistency  m  the  selection  of  key  parameters 
when  complex  models  are  used  [Beck,  1987;  Feddersen  et 
at.,  2004].  Water  depth  is  one  key  parameter  in  numerical 
wave  models  for  shallow  water  studies  and  a  second  one  is 
the  bottom  roughness,  which  contributes  to  both  bottom 
friction  and  dissipation. 

[?|  This  study  focuses  on  retrieving  the  bottom  roughness 
length  kfy  using  an  inversion  method  based  on  a  numerical 
wave  model  Determining  the  appropriate  distribution  of 
bottom  roughness  for  a  gi  ven  wave  field  w  ith  a  numerical 
wave  model  requires  the  following  conditions  to  be  met: 
( i )  the  model  equations  must  capture  the  physical  relation¬ 
ship  between  bottom  friction  and  the  wave  Held:  (2}  it  is 
necessary  to  have  a  dense  set  of  observations;  and  (3)  the 
solution  must  be  unique  for  the  selected  key  parameter.  The 
numerical  wave  model  must  include  a  bottom  friction 
dissipation  term  in  the  model  equations  that  explicitly 
represents  (he  nonlinear  damping  of  wave  energy  caused 
by  bottom  roughness  (condition  1),  We  use  a  numerical 
wave  model  to  produce  a  database  that  represents  a  dense 
set  of  observations  (condition  2).  We  will  meet  condition  3 
by  demonstrating  that  the  inversion  works  for  small  bathy¬ 
metric  errors  because  the  influence  of  water  depth  on  the 
wavefteid  is  distinct  from  that  of  bottom  friction. 

f(»]  This  paper  describes  an  inversion  method  for  deter¬ 
mining  the  bottom  friction  Held  in  a  coastal  area  for  which 
the  bathymetry  is  known,  The  purpose  oJ  the  technique  is  to 
learn  more  about  the  bottom  using  the  physics  of  energy 
dissipation  by  bottom  friction.  The  method  utilizes  a 
numerical  wave  model  in  an  iterative  procedure  similar  to 
previous  studies  |e.g..  Kennedy  et  ai ,  2000;  Narayanan  et 
at.v  2004).  The  Kennedy  et  aL  [2000]  approach  uses 
observations  of  wave  speed  to  make  local  corrections  to 
the  bathymetry.  In  the  present  study*  I  he  global  impact  of  a 
modified  variable  on  the  observations  is  contained  within  an 
influence  matrix,  which  permits  nonlocal  effects  to  be 
included  in  the  inversion.  We  will  demonstrate  that  this 
approach  can  be  used  if  the  key  parameters  are  well  defined. 
This  study  will  focus  on  wave  height  because  it  can  be 
measured  by  both  in  situ  and  remote  observation  methods. 
This  inversion  method  is  intended  for  use  with  dense  data 
sets  from  remote  sensing,  including  satellites  and  coastal 
radar  systems  |c,g„  Haas,  2007).  However,  no  such  obser¬ 
vations  were  available  for  this  study*  so  these  data  sets  will 
be  represented  by  a  series  of  control  runs  using  a  numerical 
wave  model.  The  inversion  procedure  wilt  then  attempt  to 
reproduce  these  wave  Helds  while  recovering  the  original 
bottom  roughness  distributions. 

[?]  Section  2  describes  the  numerical  wave  model  used  in 
this  study  and  the  model  setup  for  the  different  eases 
examined.  Section  3  presents  the  inversion  method  and 
section  4  discusses  the  results  for  retrieved  bottom  friction 
fields  Section  5  examines  the  sensitivity  of  the  inversion  to 
noise  and  unknown  errors  in  bathymetry  and  section  6 
discusses  some  of  the  issues  inherent  in  applying  this 
approach. 

2,  Wave  Simula  (ions 

[s]  The  approach  described  in  this  paper  docs  not  depend 
on  the  wave  model  used  but  it  is  necessary  that  the  model 


contains  the  key  parameter  of  interest  We  use  the  SWAN 
(Simulating  Waves  Nearshore)  model,  which  includes  terms 
for  refraction,  reflection  caused  by  currents  and  water  depth, 
wave  breaking,  dissipation  by  bottom  friction  and  white- 
capping.  wave- wave  interaction,  and  local  wind  generation. 
Recent  versions  of  SWAN  also  calculate  diffraction,  but 
they  were  not  used  in  this  study.  SWAN  incorporates  depth- 
induced  breaking  but  it  is  turned  off  for  this  study.  White 
capping*  wind  input*  and  nonlinear  interactions  are  also 
turned  off. 

[v]  The  SWAN  wave  model  represents  surface  waves 
with  the  two-dimensional  wave  action  density  spectrum 
where  o  is  the  relative  frequency  and  0  is  the  wave 
direction  [Baoij  et  at .*  1999|.  The  wave  spectrum  is 
described  by  the  spectral  action  balance  equation: 


where:  V  the  gradient  operator  in  x,  y,  th  and  a:  C  the 
group  velocity;  and  the  source  term  5'  is  given  by  S  S,>f  + 
Satfnrn  For  a  discussion  of  the  other  source  terms  Suf/t^  the 
interested  reader  is  referred  to  Bootj  et  aL  [ )  999).  The  bottom 
friction  term  S!ff  is  defined  by: 


Sm  (a,  B)  s=  — Q 


g’sinh (kd) 


(2) 


where:  g  -  the  gravity  constant:  A  =  wave  number:  d  -  water 
depth;  and  En  -  the  energy  density.  The  drag  coefficient  C(tf 
is  given  by  Chf  fw  -7^  Urtm  where  is  a  bottom  orbital 
velocity.  The  wave  friction  factor/,  is  found  by  solving  the 
following  equation: 


4\/7^ 


togt. 


1 

4 


-  0 .08  1  logl0 


kfy 


13) 


[in]  Equation  (3)  contains  the  near-bottom  excursion 
amplitude  ati  and  the  bottom  roughness  kh.  The  bottom 
friction  formulation  used  in  this  study  is  the  eddy  viscosity 
model  of  Madsen  et  aL  [1988].  Three  values  of  the 
roughness  length  kit  are  used:  0.005  m*  0.05  m,  and 
0.07  m.  Sensitivity  tests  with  the  wave  climate  used  herein 
showed  that  for  larger  values  of  A*,  there  is  an  effective  cap 
on  the  bottom  friction  term  for  the  Madsen  formula.  A 
roughness  length  kt,  of  0.05  m  is  used  as  the  background 
value  in  all  of  the  simulations  because  ii  is  a  typical  value 
associated  with  ripples  on  the  inner  shelf 
[11]  SWAN  integrates  (l)  with  finite  difference  schemes 
for  time,  geographic  space*  and  spectral  space.  A  cun  il inear 
grid  (Figure  I )  is  used  for  the  prognostic  wave  calculations 
by  SWAN  and  the  wave  He  Ids  are  output  on  a  uniformly 
spaced  grid  (hereinafter  termed  output  grid  or  OG)  using 
bilinear  interpolation,  which  is  incorporated  within  SWAN, 
The  OG  is  representative  of  remotely  sensed  or  measured 
wave  fields.  Nodes  on  this  grid  are  denoted  i„  and  /,  along 
the  .x  (easting)  and  v  (northing)  axes*  respectively.  A 
uniform  output  grid  is  more  convenient  for  processing  the 
wave  height  and  wave  height  error  fields.  The  use  of 
different  calculation  and  output  grids  introduces  some 
interpolation  errors,  which  are  considered  acceptable  for 
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the  simple  vvavelldds  considered  herein.  The  wave  height 
error  analysis  and  model  correction  compulations  are 
performed  on  a  curvilinear  grid  (hereinafter  called  the 
analysis  grid  or  AG).  The  analysis  grid  must  be  identical 
to  or  comprise  a  subset  of  the  computational  grid.  In  ibis 
study,  the  curvilinear  analysis  grid  is  created  by  using  every 
third  row  and  column  from  the  curvilinear  computational 
grid:  thus  the  number  of  cells  is  decreased  by  u  factor  of 
nine,  which  does  not  introduce  interpolation  error  but  there 
is  a  loss  of  information  on  the  AG.  Nodes  in  this  grid  are 
denoted  by  ia  and  ju  along  the  easting  and  northing  axes, 
respectively, 

[n]  The  setup  for  SWAN  is  the  same  lor  all  experiments* 
except  for  the  bottom  roughness  distribution  and  bathyme¬ 
try.  The  computational  grid  spans  88  km  along  the  x  axis 
with  1 7b  cells  having  a  uniform  spacing  of  500  m.  and 
25  km  perpendicular  to  the  coast  (r  axis)  with  45  cells 
whose  dimensions  decrease  from  1000  m  at  the  seaward 
margin  to  100  m  at  the  coast  (Figure  l).  The  resulting 
bottom  gradient  of  1.5  x  10  1  is  between  that  of  the  U,  S. 
middle- All  antic  shelf  (less  than  1.1  x  3  0  )  and  that  of  the 

Gulf  of  Mexico  (approximately  6  x  111  4).  The  eo input a- 
tional  grid  has  a  stretched  v  coordinate  because  wave 
damping  is  insignificant  in  deeper  water  and  thus  high 
resolution  is  unnecessary  there.  A  uniform  grid  with  a  cell 
size  of  200  ni  is  used  for  output  from  SWAN,  The  open 
boundary  forcing  consists  of  a  JON  SWAP  spectrum  applied 
uniformly  on  the  southern  boundary  with  a  peak  enhance¬ 
ment  factor  of  3.3.  The  mean  wave  lieiglu  is  I  m  and  the 
wave  source  direction  is  oriented  30*  from  shore  normal 
(coining  from  south  southeast).  There  is  no  smooth  and 
realistic  way  to  specify  a  lateral  boundary  condition  for 
Lhis  grid  orientation  in  SWAN.  The  western  and  eastern 
boundaries  use  no  wave  forcing,  which  means  that  the 
predicted  results  are  invalid  for  x  greater  than  60  km  (see 
Figure  I  ).  A  directional  resolution  of  5°  is  used.  Depth- 
induced  breaking  is  represented  in  the  wave  model  and 
would  not  be  a  problem  if  it  were  included  in  both  the 
control  and  inversion  runs;  however,  it  is  i umed  o IT  for  this 
study  because  the  spatial  resolution  of  the  grid  (e,g..  grid 


depths  of  0,4  m,  1.4  m,  and  2.6  m  within  the  surf  zone)  is 
inadequate.  The  model  is  used  to  simulate  swell  conditions 
only  (no  wind  input,  w  hilecapping  or  nonlinear  interactions) 
and  thus  the  compulations  are  stationary'.  We  chose  to  use 
three  iterations  to  compute  the  stationary  wave  Held  using 
the  implicit  scheme  used  by  SWAN. 

[n]  This  study  focuses  on  a  bottom  friction-dominated 
problem,  which  is  reflected  in  the  included  model  physics. 
This  simplified  model  reduces  the  numerical  problem  to  one 
with  only  one  dominant  physical  process,  wave  damping  hv 
bottom  friction.  Thus  the  key  parameter  is  well  known  and 
problems  of  parameter  identification  arc  minimized.  Both 
refraction  and  shoaling  are  included  in  the  model  and  are 
expected  to  be  represented  accurately.  In  order  to  increase 
the  signal-to-noise  ratio  in  the  experiments,  winds  and  other 
coastal  forcing  fields  are  not  included  because  they  would 
increase  environmental  noise  that  is  not  explicitly  addressed 
in  this  study.  A  background  (i.e.,  regional)  bottom  rough¬ 
ness  is  assumed  known  in  order  to  increase  the  speed  of 
convergence  to  solution,  but  it  is  not  required  by  the 
method. 

[  1 4  ]  The  re  are  t  wo  l  y  pes  o  f  geo  log  i  ea  I  feat  u  re  s  u  sed  i  n  1  he 
simulations,  bottom  type  (friction  features)  and  bathymetry. 
The  friction  features,  which  are  on  the  order  of  100  km2  in 
area,  represent  variations  in  bottom  characteristics  such  as 
bed  forms  and  sediment  type.  The  spatial  scale  of  sediment 
heterogeneity  varies  bul  mud  banks  can  be  as  much  as 
1800  km2  in  extent  [Anthony  et  ai,  2002]  whereas  individual 
ripples  are  much  smaller.  The  value  of  kh  used  for  these 
features  in  the  model  represents  spatial  averages  of  ripple 
Helds.  The  friction  features  are  simplifications  of  these 
bottom  types  that  use  only  binary  friction  values  (e.g,,  kh  - 
0.07  m  or  0.005  m).  This  simple  structure  is  unknown  to  the 
inverse  method,  which  does  not  utilize  information  about  the 
parameter  distribution  in  its  solution.  Therefore  the  inverse 
solutions  comprise  a  range  of  values  rather  than  binary.  The 
bathymetric  depressions  represent  features  like  drowned 
river  valleys  and  the  elevated  areas,  which  have  horizontal 
dimensions  of  less  than  3000  m,  represent  shoals  and  sand 
ridges  [McBride  and  Mostmw  1 00 1 J 

[15]  A  total  of  eight  eases  are  discussed  in  this  report 
(  fable  I ).  These  examples  illustrate  how  well  the  inversion 
can  recover  an  unknown  bottom  roughness  distribution. 
Case  (l  )  is  a  trivial  simulation  that  has  no  variations  in 
bottom  friction  or  the  bathymetry  in  Figure  1.  It  is  used  to 
compute  the  influence  matrix  as  discussed  in  section  3.1. 
Cases  (2)  through  (5)  examine  different  kh  distributions 
(Figure  2)  that  represent  patches  of  cither  rougher  {kfl 
(M)7  m)  or  smoother  (kb  0.005  m)  bottom  on  an  otherw  ise 
uniform  seafloor  (kb  ^0.05  m).  Cases  (5)  and  (6)  introduce 
random  noise  into  the  wave  fie  Id  to  determine  if  the  inver¬ 
sion  can  recover  a  bottom  roughness  field  ilial  includes  a 
low-friction  patch.  Case  (7)  incorporates  an  unknown  em¬ 
bay  men t  that  represents  a  large  error  in  bathymetry*  and 
case  (8)  includes  unknown  localized  bathymetric  highs 
(amplitude  =  0.5  2  m)  as  well  as  a  patch  with  reduced 
bottom  roughness, 

3*  Inversion  Method 

ftf/]  The  inversion  described  in  this  paper  relies  on 
computing  a  large  number  of  solutions  using  a  numerical 
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Table  t,  Boundary  Conditions  and  Parameters  Used  for  SWAN  Runs 


l  .IvC 

Shelf  Gradient  and 

Bathymetry 

kh*  m 

Wave  field'1 

r 

turns 

(MIS 

OBC:  It* 

1  m  from  SSf- 

2 

0.0015 

Background  0.05  ni  and  feu  lure  11  005  m 

OBC:  tl» 

1  m  from  SSE 

y 

0,00 1 5 

Background  0,05  m  and  feature  0  07  m 

OBC:  Hs 

1  m  from  SSE 

4 

0,0015 

Baekgmund  0.05  m  and  IWO  features  0.07  in 

OBC:  liN 

1  m  from  SSf 

5 

0.0015 

Background  0.05  m  and  feature  0.005  m 

Case  2  field 

with  random  error  of  ±11.03  in 

6 

0.1X1 1 5 

Background  0.05  ni  and  feature  0  005  m 

Case  2  field  with  random  emir  of  40. 1 2  m 

7 

Background  0,00 1 5  and  embay  men! 

0,05  m 

OBC:  Ik 

1  m  from  SSC 

S 

Bnckfiround  0,00 1 5  with  shoals 

Background  0.05  m  and  feature  0,07  m 

OBC:  Ik 

1  m  from  SSF 

*OHC  is  open  boundary  cmidiiinn,  .md  SNE  is  south  southeast. 
Yasc  I  is  only  used  for  calculating  IM 


wave  model  and  II riding  the  optimal  solution  by  minimizing 
the  error  in  the  model-predicted  significant  wave  height. 
The  penally  in  terms  of  slower  calculations  is  considered 
acceptable  because  it  is  possible  to  improve  the  numerical 
technique  later  to  increase  efficiency.  One  advantage  is 
that  the  procedure  is  independent  of  the  wave  model  used. 
A  similar  approach  was  used  by  Kennedy  et  at.  [2000] 
lor  determining  inverse  depths.  The  algorithm  for  removing 
the  roughness  field  consists  of  the  following  sLcps. 

1 1  ?]  1.  Collect  dense  measurements  of  the  wave  height 
field  over  the  area  of  interest  at  one  time. 

[ih]  2-  Estimate  the  background  bottom  roughness  field 
for  ihc  region  of  interest;  convergence  is  faster  for  an 
accurate  estimate, 

[wj  3.  Calculate  the  normalized  wave  height  change 
relative  to  Case  1  (the  initial,  homogeneous  estimate  of 
the  friction  field)  on  the  output  grid  associated  with  a  small 
perturbation  of  the  bottom  friction  at  each  point  on  the 
analysis  grid.  This  matrix,  which  has  a  number  of  columns 
and  row  s  equal  to  the  number  of  cells  on  the  OG  and  AG* 
respectively,  is  the  Influence  Matrix,  IM. 

[10]  4.  Run  the  numerical  wave  model  with  the  initial 
bottom  friction  field  (or  field  updated  in  step  5)  to  find  the 
stationary  wavcfield, 

t  j  5,  Compare  the  measured  and  predicted  wine  heights 
and  use  this  result  along  with  IM  to  update  the  bottom 
friction  field. 

[zz]  6.  If  the  maximum  wave  height  error  is  above 
a  specified  tolerance,  go  to  step  4:  otherwise,  end  the 
sequence. 

[m]  This  algorithm  contains  two  key  procedures,  steps  3 
and  5  The  Influence  Matrix  computed  in  step  3  includes  the 
impact  of  a  change  in  bottom  friction  at  a  given  point  on 
every  other  point  wiihin  the  model  grid.  It  is  computed  only 
once  using  the  initial  estimate  as  the  friction  field.  The 
friction  field  is  updated  in  step  5  by  using  the  column  in 
IM  for  the  cell  on  the  OG  with  the  largest  error.  Steps  4 
through  6  comprise  the  analysis  cycle.  The  number  of 
analysis  cycles  required  is  different  for  each  of  the  cases 
discussed  in  this  paper, 

3,1,  Influence  Matrix  Calculation 

[m]  Tlie  Influence  Matrix,  which  is  calculated  in  step  3 
of  the  inversion  algorithm,  contains  no  x  mo  columns  and 
na  x  ma  rows*  where  no  and  mo  are  the  number  of  cells 
along  the  x  and  y  axes  on  the  output  grid  and  na  and  mu  are 
the  axes  dimensions  ofihe  analysis  grid.  Each  point  (/,,,  /„) 
on  the  OG  corresponds  to  a  column  n  in  IM.  where  n 
1. ,  .no  x  mo,  and  every  point  on  the  AG  corresponds 


to  a  row  m  in  IM,  where  m  1. .  na  x  mu.  The  elements  of 
IM  are  then  found  by  repeating  the  following  steps  for  each 
point  on  the  AG,  or  row  m  in  IM. 

(z«]  l.  The  bottom  roughness  length  at  cell  (/,„  jit) 
corresponding  to  row  m  is  modified  from  the  initial  estimate 
by  multiplying  by  a  factor,  -\kh.  The  value  of  y  should  be 
sufficient  to  have  a  noticeable  effect  (or  influence)  on  the 
wave  height  field.  We  use  -7  =  0.1  on  the  basis  of  a  number 
of  test  cases. 

[zit]  2  A  new  wave  field  is  computed  with  the  numerical 
wave  model  using  the  modified  friction  field. 

[27]  3.  The  normalized  change  in  the  predicted  wave 
height  Hrn  on  the  outpui  grid  (or  column  n  in  IM)  caused 
by  the  modified  kh  at  cell  U‘„v/rJ)  is  then  given  by 

.»  ti\ M 

hmM  - - 77— -  “  I  .no  x  mo  (4) 

©Max 


w  here:  H  (  „  the  wave  height  predicted  by  Case  I  ( I  igure  3a ) 
for  each  cell  on  the  OG  and  E\ |A\  maximum  value  of  E, 
The  normalized  change  ranges  from  0  to  1.  Note  that  E  is 
positive  because  it  represents  the  strength  of  the  influence 
only. 

[zh]  After  completing  these  three  steps  for  all  rows  {m 
I . .  jta  x  /Ha),  which  correspond  to  the  points  on  the 
analysis  grid,  all  of  the  elements  in  IM  have  been  calculated: 


fc'u 

/■;,  2 

fe'u 

■  ■  t+l 

&I,I 

Ell 

Bii  ■ 

■  1 

Eu 

fc'jj 

*-  EyjHtvmtt 

t'*Ul  '  flrf.l 

&rt.t  \ m.i 

2  K,u -  wm  ,  \  ■ 

‘  ■  E,U,  „  W|)||M  pm,  I 

(5i 


[29]  The  waveficlds  computed  by  the  numerical  model 
are  dependent  on  the  boundary  conditions  and  bathymetry; 
therefore  it  is  unnecessary  to  recompute  IM  unless  these  are 
intentionally  changed.  For  application  to  the  coastal  ocean, 
IM  would  be  calculated  for  a  steady  wavcfield  using 
available  bathymetry  but  it  would  need  to  be  recomputed 
if  the  wavcfield  changed  during  the  sludy  interval.  The 
numerical  wave  mode!  is  run  na  x  ma  times  (H26  for  the 
AG  used  in  this  study);  of  course*  these  are  stationary 
wav  efields,  which  greatly  reduces  the  computations  required. 

[jo]  Each  row  of  JIM  contains  the  normalized  wave 
height  changes  on  the  output  grid  from  (4)  for  the  specified 
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Figure  2.  The  (riel ion  coefficients  kft  used  for  wave 
simulations,  (a)  Case  2,  feature  kfi  U.005  m,  (b)  Case  3, 
feature  k{)  0.07  mt  and  (c)  Case  4,  both  features  have  kh  = 

0.005  ni.  The  background  kh  0.05  m.  The  coastline  is  on 
the  north  side  of  the  grid  (top  of  page). 


bottom  friction  change  at  a  single  point  on  the  analysis  grid. 
The  choice  of  I  lie  analysis  and  output  grids  is  not  important 
to  the  inversion  method  itself.  The  grid  attributes  are 
selected  for  convenience  of  computation  and  the  desired 
accuracy.  We  have  retained  the  multiple  grid  nomenclature 
for  generality.  The  ordering  of  grid  points  in  IM  is  not 
important  because  the  normalized  changes  E  arc  mutually 
independent. 

3.2.  Updating  the  Friction  Field 

[a  i  J  After  a  new  wave  field  has  been  computed  in  step  4 
of  the  analysis  cycle,  the  wave  heights  are  compared  to  the 
observations  (control  runs  in  these  examples)  in  order  to 
modify  the  bottom  roughness  length.  This  step  (step  5  of  the 
inversion)  is  accomplished  using  the  following  procedure, 
which  is  repealed  until  the  maximum  wave  height  error  is 
within  a  specified  tolerance. 

[jtt]  1 .  The  wave  height  error  at  each  point  on  the  OG  is 
given  by 

Eno=HAc-H<  (6) 

where:  H ,  is  the  significant  wave  height  from  the  control 
run  (c  -  2  through  8;  see  Table  1);  and  H A%  is  the  wave 
height  computed  by  the  wave  model  with  the  retrieved  kft 
field.  The  initial  estimate  of  A*  is  used  on  the  first  iteration. 
Note  that  positive  values  indicate  overprediction.  The  sign 
of  the  error  is  mapped  to  the  AG  and  saved  as  yt  for  use  in 
step  3  below, 

[m]  2,  Noise  is  removed  from  the  Efio  field  with  a  3  x  3 
(2-D  nine-point)  filler,  and  the  cell  with  the  maximum  wave 
height  error  magnitude  (/,*.,  jtHX  which  corresponds  to 
column  ne  in  IM,  is  found, 

[i4|  3.  Noting  that  every  point  (iUi  ii{)  on  AG  corresponds 
to  a  row  m  in  IM,  the  roughness  length  for  all  grid  points  on 
the  AG  is  updated  from  the  previous  analysis  cycle  by 

kL  -  kKm  +  AkL  m  -  1  ■ ««  *  m“  (7) 

where:  p  -  iteration  or  analysis  cycle:  and  „p  change 
in  roughness  length,  which  is  given  by 

Ak^  —  a  tit  =  I  ,tm  x  ma  (S) 

where  a  is  a  constant  with  the  same  units  (meters)  as  the 
key  parameter  being  retrieved.  A  larger  value  of  a  means 
more  modification  to  the  friction  field  with  each  analysis 
cycle.  The  sign  of  the  change  is  obtained  from  y/ 

[?>]  Here  we  assume  a  linear  relationship  between  the 
changes  in  bottom  friction  and  wave  height  F aster  conver¬ 
gence  might  be  achieved  through  the  use  of  a  more  accurate 
nonlinear  relation,  particularly  when  coupled  to  a  relaxation 
scheme.  The  constant  a  was  evaluated  in  a  sensitivity  study 
on  Case  4  with  three  param  clerical  ions:  (a)  a  constant  value 
of  0.01  m;  (b)  using  a  decreasing  step  function;  and 
(c)  decreasing  magnitude  with  cycle  number.  The  maximum 
error  Eff0  was  reduced  to  0.05  m  in  (>80  cycles  for  a 
constant  value  and  in  fewer  than  500  cycles  for  (c).  which 
had  the  fastest  convergence,  A  constant  value  of  0.01  m  is 
used  in  this  paper  in  order  to  avoid  changing  the  parame¬ 
terization  for  different  cases.  The  most  important  factor  in 
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FiRure  3.  The  wave  height  from  the  control  runs  (synthetic  observations),  (a)  Case  I,  homogeneous 
bottom  friction  usevl  only  to  cal  at  late  1M,  (b)  Case  2,  a  single  feature  with  decreased  bottom  friction, 
(c)  Case  3,  A  single  feature  with  increased  bottom  friction*  and  (d)  Case  4,  two  features  with  decreased 
bottom  friction.  The  coastline  is  on  the  north  side  of  the  grid  (top  of  page). 


finding  the  change  in  bottom  roughness  from  18)  is  the 
value  of  elements  i/u,  ttt.)  in  IM  for  (m  L  na  x  ma).  This 
factor  reflects  die  physical  influence  of  bottom  friction,  us 
represented  by  the  wave  model  at  cell  K  on  all  other 

cells  within  the  physical  domain, 

4*  Inversion  Results 

[>]  This  section  presents  the  results  of  applying  the 
inversion  method  described  above  to  the  cases  listed  in 
Table  J.  We  will  refer  to  simulations  using  known  friction 
fields  as  control  inns,  which  serve  as  proxies  for  the  dense 
data  sets  that  would  ideally  be  available  from  remote 
sensing  We  present  results  for  Cases  2  through  H  because 
Case  I  is  only  used  for  computing  IM.  The  convergence 
behavior  of  the  algorithm  is  evaluated  by  comparing  the 
significant  wave  heights  computed  during  the  analysis  to 
those  from  the  control  run  (synthetic  observations).  The 
predicted  wave  heights  are  very  sensitive  to  errors  in  bottom 
roughness  and  even  small  errors  can  indicate  problems, 
requiring  examination  of  the  retrieved  roughness  fields.  The 
accuracy  of  the  solutions  will  be  examined  by  comparing 


the  retrieved  fields  to  the  known  fields  from  the  control 
runs. 

4, 1,  Waves  Over  an  Elliptical  Low-Friction  Area 

[37]  The  bottom  roughness  field  for  Case  2  (Figure  2a) 
consists  of  a  uniform  friction  field  w  ith  kh  =  0.05  m,  which 
represents  ripples  that  are  0.05  m  high*  and  an  elliptical 
feature  with  reduced  bottom  dissipation  {kh  0.005  m)  that 
extends  from  a  water  depth  of  12  m  to  a  depth  of  20  rn  In 
the  ocean  a  feature  like  ihis  could  be  due  to  finer  sediments, 
which  would  have  smaller  ripples.  The  inversion  for  this 
case  is  presented  in  detail 

[isj  The  control  run  wave  heights  (Figure  3b)  decrease 
landward  away  from  the  feature  because  the  increase  in  Us 
by  shoaling  is  negated  by  bottom  friction.  This  cross-shore 
distribution  is  consistent  with  observations  from  a  natural 
beach  with  similar  incident  waves  [Efgar  et  ai  *  1997].  The 
decreased  dissipation  over  the  feature  results  in  higher 
waves  closer  to  the  coast. 

[w]  A  plot  of  Ejio  (equation  (6))  a  Her  one  analysis  cycle 
(Figure  4a)  shows  that  the  maximum  occurs  at  the  shallow 
end  of  the  feature.  The  errors  are  negative  over  the  feature 
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itself  because  of  its  much  lower  roughness  <0.005  m)  mid 
resulting  higher  waves  in  the  control  run,  The  values  of  E 
found  from  (4)  for  column  nv  of  IM  (Figure  4b)  indicate 
where  a  change  in  bottom  friction  will  have  the  greatest 
influence  on  the  wave  height  at  For  example.  E  is 

near  zero  far  from  the  friction  feature  and  in  deep  water,  and 
large  at  the  shallow  end  of  the  feature.  The  largest  values  of 
E  are  oriented  south -southeast  because  bottom  dissipation 
impacts  the  wave  field  on  the  leeward  si  tie  of  a  friction 
feature.  The  normalized  wave  height  changes  in  Figure  4h 
are  used  in  (8)  to  calculate  the  adjustment  of  the  friction 
field  A kh  after  the  first  analysis  cycle  (Figure  4c).  The 
largest  modifications  occur  over  the  center  of  the  feature. 

[■to]  The  maximum  Efio  (Figure  5)  decreases  to  less  than 
0.01  m  (approximately  1%  for  the  l  m  waves  used  in  this 
study)  after  fewer  Ilian  50  analysis  cycles  but  the  inversion 
continues  until  the  error  is  below  a  specified  tolerance.  The 
im enbon  stops  after  261  cycles  with  a  calculated  significant 
wave  height  field  (Figure  6a)  that  is  indistinguishable  from 
the  control  run  (Figure  3b)  and  the  values  of  Et/0  arc 
near  zero  (Figure  6b).  indicating  that  the  inversion  has 
converged, 

[4i ]  The  recovered  kh  field  (Figure  6c)  is  dose  to  the 
original  (Figure  2a),  but  the  values  range  from  a  minimum 
of  0.005  m  to  the  background  value  (0.05  ml  whereas  the 
control  run  values  have  a  binary-  distribution  (either  0,005  m 
or  0.05  m),  The  normalized  mean  error  between  these 
values  and  the  original  field  is  —4%  (kly  overpredicted) 
and  the  standard  deviation  is  45%.  Note  that  the  retrieved 
bottom  roughness  is  less  accurate  in  deeper  water 

4.2.  Waves  Over  an  Elliptical  High-Friction  Area 

[4zJ  The  friction  feature  lor  Case  3  (Figure  2b)  is  similar 
in  plan  to  Case  2  hut  kft  is  increased  to  0.07  m.  A  larger 
value  of  kft  cannot  be  used  because  the  friction  term  in 
SWAN  effectiv  ely  limits  the  bottom  roughness  length.  The 
effect  of  this  feature  on  the  wave  height  (Figure  3c)  is  weak 
but  noticeable. 

[  n]  The  maximum  wave  height  error  (Figure  7)  decreases 
from  approximately  0.02  m  after  one  cycle  to  less  than 
0.0006  m  after  41  cycles.  The  error  is  positive  because  the 
initial  estimate  for  kft  (0.05  m)  causes  overprediction  of  the 
wave  heights  over  the  patch  The  dependence  of  the  maxi¬ 
mum  E//0  on  cycle  number  is  similar  to  Case  2  (Figure  5 )  but 
the  solution  converges  much  faster  because  of  the  reduced 
effect  of  A*  0.07  m.  This  result  can  be  seen  by  comparing 
the  analysis  wave  heights  (Figure  8a)  to  the  control  run 
(Figure  3c), 

[44]  The  final  bottom  roughness  field  (Figure  8b)  is 
similar  lo  the  initial  field  shown  in  Figure  2b.  The  mean 
error  for  the  retrieved  kh  field  is  less  than  -1%  and  the 
standard  deviation  of  the  error  is  less  than  4.8%.  This  result 
is  consistent  with  the  small  errors  in  the  predicted  wave 
heights. 

43.  Waves  Over  Two  Cross-Shore  Low-Friction  Areas 

[45 1  There  arc  two  low- friction  (kh  0.005  m).  shore 
normal  features  in  Case  4  (Figure  2c).  This  pattern  is  more 
representative  of  a  typical  inner  shelf  and  sho  reface  than  the 
previous  cases  but  ii  also  increases  the  difficulty  of  the 
solution.  This  bottom  toughness  field  has  the  greatest  c  fleet 
on  the  control  run  wave  heights  (Figure  3d)  of  the  three 
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Figure  5,  Case  2;  the  maximum  value  of  the  filtered  wave 
height  error  bl(j(h  us  a  function  of  analysis  cycle. 


cases  discussed  so  far,  Wave  heights  increase  at  two  points 
along  the  shore  because  of  the  features:  however,  neither  the 
features  nor  the  wave  lie  Id  are  symmetric. 

fi*J  The  maximum  En()  a  Her  one  cycle  is  0.3  m  (  figure  9), 
which  is  much  greater  than  for  the  previous  cases.  The 
maximum  occurs  within  the  western  feature.  The  values  of 
£  from  IM  corresponding  to  cell  (/,„.,  jiu  )  (not  shown) 
indicate  that  correlation  distances  are  much  shorter  for  these 
shoretine-uuaehed  features  than  for  offshore  features.  This 
is  consistent  with  the  physical  problem;  in  shallower  water, 
the  ratio  between  near- Held  wave  energy  dissipation  and 
far- field  dissipation  becomes  larger:  thus  the  changes  in 
bottom  friction  that  are  used  in  computing  l\1  do  not  have 
a  far-reaching  influence. 

[41  j  The  maximum  Euo  decreases  to  less  titan  0.03  m 
after  800  cycles  but  the  inversion  continues  until  cycle 
1000.  The  largest  errors  (Figure  10a)  are  in  shallow  water 
east  of  lire  features.  Positive  values  of  EJf0  result  from 
overcome  lion  of  kh  where  no  differences  in  bottom  friction 
actually  existed.  This  is  caused  by  the  choice  of  the  AG, 
which  spreads  the  correction  from  IM  over  more  cells  on 
the  OG,  The  large  error  (^3%)  suggests  that  the  retrieved  kh 
field  should  he  examined  in  detail. 

[4*|  The  retrieved  friction  field  (Figure  1 0b)  resolves  the 
two  features  but  their  shapes  have  not  been  fully  recovered 
offshore  because  the  values  of  E  contained  in  IM  are 
relatively  small  in  deeper  water.  The  mean  normalized  error 
for  kft  is  59%  and  the  standard  deviation  is  183%.  This 
result  suggests  that  even  small  errors  in  wave  height  can 
indicate  a  problem.  Furthermore,  the  mean  error  is  biased  by 
the  solution  for  waier  depths  greater  than  15  m,  where 
bottom  friction  has  minimal  effect  on  the  incident  waves. 
This  problem  demonstrates  the  limitations  of  using 
measured  wave  heights  to  infer  bottom  friction  coefficients 
in  deep  water  where  wave  heights  are  insensitive  to  bottom 
friction.  The  maximum  wave  height  error  does  not  decrease 
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Figure  ft.  Case  2:  results  from  the  inversion  analysis  at 
cycle  261 .  (a)  The  analysis  wave  height  field,  (b)  plot  of  the 
filtered  wave  height  error  EfHh  and  (c)  contour  plot  of  kh 
values.  The  coastline  is  on  the  north  side  of  the  grid  (lop  of 
page). 
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Figure  7.  Cm1  3:  the  maximum  value  of  ihe  filtered  wave 
height  emir  Eitth  as  a  funetion  of  analysts  cycle. 


below  0.02  m.  which  is  an  order  of  magnitude  greater  than 
for  the  previous  runs.  It  is  thus  important  to  note  that  values 
of  Euo  greater  than  1%  are  an  indication  of  poor  conver¬ 
gence  of  the  inversion,  and  an  accompanying  uncertainty 
in  ihe  retrieved  bottom  friction  Held,  Nevertheless,  the 
recovered  kft  field  (Figure  1 0b)  is  reasonable  tor  depths 
shallower  than  15  m. 

5.  Inversion  Sensitivity  and  Uniqueness 

[^\  The  tests  discussed  in  this  section  use  the  same 
boundary  condition  as  in  Cases  I  through  6  and  assume 
the  same  bathymetry';  therefore  IM  does  not  need  to  be 
recalculated*  We  will  thus  evaluate  the  impact  of  uncertainty 
in  the  measured  wave  field  (Cases  5  and  6)  and  examine  the 
inversion  result  for  bottom  friction  when  another  key 
parameter  may  be  poorly  known:  specifically,  bathymetry 
(Cases  7  and  8). 

5.1*  Measured  Wave  Height  Uncertainty 

[siij  Case  5  investigates  the  effect  of  uncertainty  in  the 
wave  measurements  by  adding  random  noise  to  the  control 
run  wavefidd.  When  random  noise  with  a  standard  devia¬ 
tion  of  0.03  m  (3%  of  wave  height  at  the  open  boundary)  is 
added  to  the  wave  heights  from  Case  2.  the  inversion 
terminates  in  53  cycles  with  a  maximum  filtered  wave 
height  error  of  less  than  0.015  in,  a  mean  error  of  ■  10% 
for  k/t,  and  a  standard  deviation  of  60%.  It  is  also  instructive 
to  find  at  what  level  of  uncertainty  the  inversion  fails.  7‘he 
inversion  stops  after  14  iterations  when  random  noise  with  a 
standard  deviation  of  12'%  is  added  to  the  wave  height  field 
(Case  6).  The  value  of  E  in  IM  is  zero  for  column  nc  and  the 
inversion  stops  because  it  cannot  reduce  the  error  at  this 
point  further  via  modification  of  the  friction  field  The  mean 
and  standard  deviation  of  the  kf>  error  for  the  entire  domain 
are  -26%  and  129%.  respectively.  These  errors  are  larger 


than  for  Case  2,  which  has  no  noise  in  the  wavefidd,  but 
they  remain  less  than  for  Case  4. 

5,2.  Bathymetry  Uncertainty 

[mJ  ft  is  important  to  examine  the  uniqueness  of  the 
inversion  method  with  respect  to  uncertainty  in  nearshore 
bathymetry  because  water  depth  is  u  key  parameter  for 
the  wave  model.  The  purpose  of  this  test  is  to  see  if  the 
uncertainly  in  water  depth  is  erroneously  retrieved  bv  the 
inversion;  that  is,  will  the  inversion  incorrectly  identify  a 
bathymetric  feature  as  a  low- friction  area? 

[s i]  The  grid  for  Case  7  incorporates  an  embayment  with 
a  mean  depth  of  5  m  (Figure  I  la)  and  has  a  uniform  bottom 
roughness  field  ikft  0,05  in).  The  inversion  assumes  that 
the  depression  is  unknown.  Such  large  bathymetry  errors  are 
unlikely  but  this  problem  will  test  the  robustness  of  the 
inversion  algorithm  for  large  errors  in  more  I  ban  one  key 
parameter.  The  predicted  wave  heights  over  the  feature  are 
similar  to  those  further  offshore  (Figure  I  lb)  because  the 
incident  waves  have  propagated  further  before  encountering 
shallow  water  where  dissipation  by  bottom  friction  is 
substantial. 
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Figure  8.  Case  3:  final  results,  (a)  The  analysis  wave 
height  field  after  4 1  cycles  and  (b)  contour  plot  of  k},  values. 
The  coastline  is  on  ihe  north  side  of  the  grid  (top  of  page). 
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Figure  9,  Case  4:  plot  of  the  maximum  value  of  the 
filtered  wave  height  emir  as  a  function  of  analysis  cycle. 


[>*]  The  inversion  terminates  a  Her  343  cycles  with  an 
error  of  almost  0  02  m  in  //$  because  ihc  maximum  wave 
height  error  occurs  at  a  location  that  is  not  sensitive 

to  changes  in  she  friction  field.  Application  of  IM  is 
therefore  not  possible  because  all  values  for  column  nt  are 
zero.  The  maximum  error  of  0,02  m  suggests  that  the 
observed  variation  from  the  control  wave  field  is  not  caused 
by  bottom  friction  only,  and  thus  the  retrieved  values  of 
should  be  examined. 

[  5j]  As  with  Case  3.  the  greatest  modification  of  the 
background  kf,  is  in  the  shallowest  water.  In  fact,  the 
retrieved  bottom  roughness  Held  (Figure  lie)  minors 
the  embay ment,  replacing  water  depth  with  smaller  values 
of  kf}.  The  errors  increase  as  l he  coast  is  approached  and  a 
minimum  kh  of  0.005  is  retrieved.  This  example  demon¬ 
strates  why  a  fundamental  assumption  of  the  inversion 
method  is  that  the  bathymetry  is  known. 


[ 55 J  We  can  further  examine  the  inversion's  sensitivity  to 
bathymetry  errors  with  a  more  realistic  example.  Case  S  uses 
a  high- friction  feature  (Figure  12a)  with  kh  0.07  m 
superimposed  on  a  uniform  background  roughness  of 
0.02  rn.  The  bathymetry  includes  multiple  shoals,  with 
amplitudes  of  0.5  2  m  in  water  depths  ranging  from  ft  to 
2ft  m,  superimposed  on  the  smooth  bathymetry  used  in 
cases  ( 1 )  through  (6)  (Figure  12b).  The  wavefield  for  the 
control  run  (Figure  13a)  has  lower  wave  heights  over 
the  friction  feature  because  of  increased  dissipation  by  bottom 
friction.  The  waves  are  higher,  however,  over  the  shoals 
because  of  shoaling  effects.  This  can  lx*  seen  by  comparison 
of  the  bathymetry  in  Figure  1 2b  and  the  wavefield  (Figure  1 3a) 
along  the  western  end  of  the  grid  at  y  =  7000  m. 

[.**]  The  predicted  wave  heights  after  31  analysis  cycles 
(  Figure  1 3b)  reflect  the  influence  of  the  friction  feature  only, 
however,  and  nol  ihc  shoals.  The  largest  values  of  E/fo 
(Figure  1 3c)  occur  in  shallow  water  over  the  feature  and  not 
where  the  bathymetry'  errors  are  located.  On  the  basis  of  our 
previous  results,  this  large  error  suggests  that  the  inversion 
is  contaminated  by  another  key  parameter  besides  bottom 
friction  (i.e.,  bathymetry). 

\$-t]  'T  he  mean  and  standard  deviation  for  the  kh  errors  are 
less  than  -1%  and  15%,  respectively,  and  comparing  the 
retrieved  bottom  roughness  field  (Figure  13d)  to  the  original 
field  (Figure  12a)  shows  that  the  largest  Ay,  error  occurs  in 
shallow  water  and  is  not  associated  with  the  unknown 
shoals.  Ibis  example  demonstrates  the  importance  of 
examining  the  distributions  rather  than  relying  on  overall 
statistical  parameters,  which  include  errors  in  deep  water  as 
well  as  shallow.  Given  this  precaution,  it  appears  that  ihc 
inversion  recovers  the  unknown  friction  feature  reasonably 
well  (tow-A>  error)  while  indicating  that  there  arc  other 
important  features  within  the  region  that  affect  wave  height 
(high  However,  as  suggested  by  case  4,  whenever  the 
wave  height  error  is  more  than  !%  il  is  necessary  to 
examine  the  retrieved  bottom  roughness  field  to  see  if  it  is 
physically  realistic. 

6,  Discussion  and  Conclusions 

[>si]  These  results  suggest  that  the  approach  works  well 
for  friction-dominated  areas  where  swell  is  the  main  source 
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(h)  predicted  values  of  significant  wave  height,  and 
(c)  contour  plot  ol '  k/Y  values.  The  coastline  is  on  die  north 
side  of  the  grid  (top  of  page). 


of  waves,  but  there  are  several  constraints  that  must  be 
considered  in  applying  it  to  realistic  problems.  Field  data 
of  sufficient  density  to  validate  the  inversion  were  not 
available  lor  use  in  this  study  and  it  was  necessary  to  use 
synthetic  observations,  which  also  reduced  the  uncertainty 
of  the  results.  As  demonstrated  by  the  results  for  random 
noise  (Cases  5  and  6),  small  measurement  errors  slow 
convergence  but  do  not  render  the  inversion  inaccurate. 
Synthetic  Aperture  Radar  (SAK)  data  can  be  processed  to 
retrieve  wave  heights  hut  errors  can  be  as  much  as  29% 
[Mmtenbroek  and  de  talk,  2000].  Shore-based  radar 
systems  [e.g,.  Gurgel  vt  tf/.,  1999]  could  also  be  useful  for 
providing  wave  measurements  analogous  to  the  synthetic 
measurements  used  in  this  study.  These  approaches  have 
problems  but  their  continued  development  suggests  that 
dense  wave  height  observations  w  ith  errors  less  than  10% 
are  not  an  unreasonable  expectation  in  the  future 

[59]  This  study  uses  a  small  computation  domain  in  order 
to  facilitate  development  and  testing  of  the  inversion  algo¬ 
rithm.  If  a  larger  grid  had  been  used  for  the  SWAN 
simulations,  it  would  have  been  difficult  to  develop  the 
algorithm  described  in  this  paper,  which  was  intended  for 
use  on  desktop  computers.  This  is  a  common  difficulty  with 
inverse  solutions  and  data  assimilation.  There  is  no  inherent 
problem  with  sealing  the  method  to  larger  or  more  complex 
prublems.  It  should  be  noted,  however,  that  the  domain  size 
used  in  this  study  is  appropriate  for  problems  of  nearshore 
dynamics  [e,g<.  Keen  et  ai ,  2003 j. 

[«i]  A  critical  element  of  the  inversion  method  discussed 
in  this  paper  is  the  Influence  Matrix,  JM.  which  reflects  the 
impact  of  wave  propagation  throughout  the  domain  as 
represented  by  the  physics  of  the  wave  model.  Thus  the 
values  of  E  display  ellipses  with  I  heir  major  axes  aligned  in 
the  direction  of  wave  propagation  for  the  friction-dominated 
dissipation  in  this  study.  For  a  bathymetry  inversion,  how¬ 
ever,  wc  would  expect  E  to  be  more  concentrated  because  of 
the  local  effect  of  water  depth  on  wave  height  through 
shoaling  Inversions  of  the  type  discussed  in  this  paper 
should  work  therefore  as  long  as  the  key  parameters 
produce  unique  Influence  Matrices.  It  is  possible  to  recal¬ 
culate  IM  during  the  inversion  procedure  in  order  to  reduce 
1  he  maximum  wave  height  error  further  hut  this  has  not  been 
investigated  yet, 

[6i ]  Equation  (H)  contains  a  dimensional  constant  o  that 
must  be  chosen  on  the  basis  of  the  magnitude  of  the  key 
parameter  being  retrieved.  A  series  of  sensitivity  analyses 
on  the  v  alue  of  a  were  completed,  and  it  was  varied  over  an 
order  of  magnitude.  The  results  of  these  tests  indicate  that 
the  magnitude  of  a  has  a  significant  effccl  on  the  speed  of 
convergence  but  only  a  minor  impact  on  the  final  solution, 
with  a  speed  up  of  28%  between  a  constant  of  0,01  m  and  a 
variable  magnitude.  A  larger  value  means  more  modi  11  ca¬ 
tion  to  tlie  friction  field  with  each  analysis  cycle,  which  is 
good  for  the  first  few  cycles.  Thus  we  can  incorporate  a  lack 
of  know  ledge  of  the  correlation  between  grid  points  into  its 
magnitude:  for  example,  if  we  are  uncertain  of  the  influence 
of  distant  effects  we  would  want  to  use  a  smaller  value. 

[02]  The  inversion  method  described  in  this  paper  is 
similar  to  traditional  data  assimilation  approaches  but  there 
are  important  differences.  For  example,  this  work  is  focused 
on  the  need  to  leant  more  about  bottom  type  rather  than 
improve  wave  forecasting  [Sherentei  and  Stone*  2003].  This 
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Figure  12*  Case  8:  friction  feature  ami  bathymetry*  (a)  Contour  plot  of  kh  values  and  (b)  bathymetry. 
The  coastline  is  on  the  north  stile  of  the  grid  (lop  of  page). 


is  a  difficult  research  subject  because  of  the  problem  of 
measuring  bottom  friction  parameters  in  the  ocean.  Inverse 
methods  that  focus  on  key  parameters  are  important  to 
identity  fundamental  relationships  between  measurable 
wave  properties  and  bottom  eharae  ten  sties*  Another  impor¬ 
tant  consideration  in  developing  an  approach  to  study  this 


problem  is  the  commitment  lo  development  of  an  adjoint 
model.  Examining  basic  eause-and-elTeei  relationships  does 
not  require  such  an  extensive  approach  and  may  even  be 
better  served  bv  occasionally  using  different  wave  models. 

[«]  This  paper  includes  a  description  of  the  inversion 
method  and  presents  the  results  of  some  simple  tests  that  are 
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intended  to  demonstrate  its  usefulness  in  friction-dominated 
inner  shelf  areas.  In  addition*  the  analysis  of  measurement 
error  indicates  that  il  is  robust  and  will  work  with  real 
observations-  However,  the  inversion  incorrectly  adjusts  the 
bottom  roughness  length  when  the  bathymetry  is  poorly 
known.  This  effect  is  less  for  random  uncertainty  than 
systematic  errors  like  ihe  embay ment  examined  in  Case  7, 
This  is  an  important  point  because  water  depth  is  another 
key  parameter  lor  numerical  wave  models.  T  his  highlights 
the  need  to  identify  the  key  parameters  in  the  model  and 
evaluate  their  sensitivity  before  applying  this  kind  of 
approach.  1  he  ability  to  isolate  bottom  friction  errors  from 
water  depth  errors  is  crucial  lo  evaluate  but  tom  sediment 
type  using  wave  observations*  but  it  may  be  possihle  to 
perform  joint  inversions  for  bottom  friction  and  bathymetry 
[Holland*  2001],  The  inversion  presented  in  this  paper  says 
nothing  about  the  unresolved  errors  caused  by  other  factors 
but  it  seems  reasonable  to  speculate  dial  a  similar  procedure 
could  he  applied  lo  bathymetric  errors  as  well.  The  useful¬ 
ness  of  the  algorithm  relies  on  ihe  treatment  of  applicable 
physical  processes  in  the  numerical  wave  model. 

[m]  This  approach  has  limitations  that  restrict  its  appli¬ 
cability  to  scientific  problems  rather  than  improving  wave 
forecasts.  If  the  predicted  wave  height  errors  are  small 
(<l%),  it  is  probable  that  bottom  friction  is  the  dominant 
process  but  in  the  presence  of  other  important  sources  of 
error  like  bathymetry,  it  is  necessary  lo  examine  the  results 
for  physical  realism.  Nevertheless,  it  can  be  a  useful 
research  tool  that,  if  used  properly,  can  help  identify  regions 
of  likely  high  dissipation  ie*g.T  cohesive  sediment)  even 
where  uncertainty  exists  in  ihe  water  depth, 
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